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Abstract. We propose an inclusive approach for calculating characteristics of secondary electrons produced 
^ , by ions/protons in tissue-like media. This approach is based on an analysis of the projectile's interaction 

^ ' with the medium on the microscopic level. It allows us to obtain the energy spectrum and abundance 

of secondary electrons as functions of the projectile kinetic energy. The physical information obtained in 
^. , this analysis is related to biological processes responsible for the irrepearable DNA damage induced by 

(-H ■ the projectile. In particular, we consider double strand breaks of DNA caused by secondary electrons and 

O ' free radicals, and local heating in the ion's track. The heating may enhance the biological effectiveness of 

electron/free radical interactions with the DNA and may even be considered as an independent mechanism 
of DNA damage. Numerical estimates are performed for the case of carbon-ion beams. The obtained dose- 
depth curves are compared with results of the MCHIT model based on the GEANT4 toolkit. 



PACS. 61.80.-x Physical radiation effects, radiation damage - 87.53.-j Effects of ionizing radiation on 
biological systems - 41.75.Aki Positive-ion beams - 34.50.Bw Energy loss and stopping power 
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L«. 1 Introduction with the medium. We hope that this approach will facil- 

iyK itate establishing a quantitative connection between the 

Q Ion and proton beams are becoming more and more pop- amount of energy deposited into the tissue and the in- 

L; ■ ular tools for cancer therapy. Ions and protons are more duced biological consequences. Another mechanism that 

. ^ [ advantageous projectiles than the now conventional pho- we discuss in this paper is the DNA damage due to local 

^^ . tons because they may cause less damage to the healthy heating of the medium. 

H : tissues surrounding tumors and thus induce fewer side ef- Regardless of the eventual mechanism of the DNA 

. .C^ , fects. More details on the history, current status and com- damage, a projectile's propagation and stopping is the 

parison of proton a nd heav y-ion therapies can be found m ^asic process in cancer therapy and therefore we start 

the review articles |l | 2 | 3 | 4j . ^^j. g^-^^jy .^^j^j^ j^g analysis. All secondary electrons and 

The method of radiation therapy consists in the in- f^^g radicals that may cause DNA damage as well as the 

activation of cancer cells via radiation impact on their medium's temperature increase are due to energy loss by 

DNA 0. The DNA may be destroyed in a number of ^^^^^ projectiles. For numerical estimates, we have taken 

ways [g. Double strand breaking (DSB) is one of the most parameters characteristic of cancer radiation therapy us- 

effective mechanisms of irreparable DNA damage. Inter- jj^„ carbon-ion beams 
actions of low-energy secondary electrons and free radicals 

with the DNA are believed to be mainly responsible for When an energetic charged projectile enters a medium, 

the DSB's [7]. We are going to build an inclusive approach it experiences a number of atomic and nuclear interac- 

yielding such microscopic characteristics as abundance of tions. Nuclear processes, mainly fragmentation reactions, 

secondary electrons and free radicals in the region as well result in the transformation of the beam particles into 

as their energy spectrum. These characteristics can be ob- new species. The nuclear processes can be reliably mod- 

tained from the analysis of the interaction of projectiles eled with the use of Monte-Carlo simulations, in particu- 

lar by means of the MCHIT model |8j9jl0j based on the 

=" E-mail: surdutov&oakland.edu; Tel: + 1-248-370-3409 GEANT4 toolkit |11I12J . Therefore in this paper, we only 

'' E-mail:solovyov@fias.uni-frankfurt.de comment on projectile fragmentation processes without 
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considering them in detail. Atomic processes, such as elas- 
tic Coulomb scattering, scattering with excitation, single 
and multiple ionization, are much more probable but less 
destructive events for the projectile, frequently occurring 
along the projectile's trajectory. Ionization of molecules 
of the target media is the main channel of the projectile's 
energy loss. This process is the most important for biolog- 
ical applications since it results in a shower of secondary 
electrons, which may inflict more DNA damage than the 
projectile itself. 

In Section [21 we study the passage of a single ion 
through the medium, analyze the impact ionization cross 
section, linear energy transfer (LET), and evaluate the 
abundance of secondary electrons and their energy spec- 
trum. Calculations are carried out for parameters typical 
for cancer therapy. They are compared with experimental 
results and the predictions of the MCHIT model for the 
depth-dose distributions. In Section [3l we estimate the 
densities of electrons and free radicals produced by the 
projectiles. Finally, we evaluate the local heating associ- 
ated with the energy loss, and estimate DNA damage due 
to thermal effects. 
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where n is the number density of water molecules (at stan- 
dard conditions n « 3.3 x lO^^cm"'^). 

In our current analysis, the angular distribution of 
emitted electrons is not important, therefore we integrate 
over the solid angle and the DDCS is replaced with a 
single-differential cross section (SDCS). Thus, the total 
impact ionization cross section, differentiated in secondary 
electron kinetic energy, da{W,T)/dW, becomes the main 
quantity in our analysis. Besides the kinetic energy of sec- 
ondary electrons and the properties of water molecules, 
the SDCS depends on the kinetic energy of the projectile 
T and its charge z. We use the semi-empirical Rudd's ex- 
pression [15] for the SDCS, which is a parametric adjust- 
ment that combines the experimental data, calculations 
within the plane wave Born approximation and other the- 
oretical models. It is given in the following form ^15j : 
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2 Ion passage through the medium 

We first describe the passage of a single ion through bi- 
ological medium. The medium is approximated by liquid 
water, which, anyhow, is the major component of most 
living tissues. The ion typically enters the medium with a 
sub-relativistic velocity (the corresponding energy range 
used in carbon-ion therapy is 100-400 MeV per nucleon) 
and then gradually slows down because of energy loss in 
inelastic processes. As mentioned above, the impact ion- 
ization is the dominant energy-loss process. Most impor- 
tantly, this process results in the production of secondary 
electrons, which are directly or indirectly responsible for 
the radiation damage of DNA. 

The inelastic cross sections are smaller at the begin- 
ning of the trajectory (at high velocities) but then increase 
and reach a maximum at the so-called Bragg peak. In the 
vicinity of this maximum, the ion looses its energy at the 
highest rate. As its velocity decreases further, the cross 
section drops because of kinematic constraints. 

In our first analysis [l3j , we have limited our consider- 
ation to non-relativistic projectile velocities. However, the 
relativistic effects are significant for 400 MeV/u ions, and 
in this work, we include relativistic corrections. As shown 
below, this results in better agreement of our calculations 
with the measured position of the Bragg peak. 

Emission of electrons (frequently called (5-electrons) in 
proton collisions with atoms and molecules has been under 
theoretical and experimental investigation for decades |14I15| 
The quantity of interest is the probability to produce N 
secondary electrons with kinetic energy W , in the interval 
dW , emitted from a segment Ax of the trajectory of a sin- 
gle ion at the depth x corresponding to the kinetic energy 
of the ion, T, within the solid angle dfi. This quantity 
is proportional to the doubly-differentiated cross section 



dW ' ^ h \h 

% 

(1 + i^if (1 + exp(a(cu, - wr'^)M)) ' 

where the sum is taken over the electron shells of the water 
molecule, ag = 0.0529 nm is the Bohr radius, R = 13.6 eV 
is the Rydberg, Ni is the shell occupancy, li is the ioniza- 
tion potential of the shell, LOi = W/Ii is the dimcnsionless 
normalized kinetic energy of the ejected electron, Vi is the 
dimcnsionless normalized projectile velocity given by 



~2r 
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where m is the mass of electron and V is the velocity 
of the projectile. In the non-relativistic case, the expres- 
sion ^ is the square root of the kinetic energy of an elec- 
tron having the same velocity as the projectile, normal- 
ized by the corresponding binding energy. The definition 
of Vi remains in the relativistic case, however, the projec- 
tile velocity V is no longer equal to \ "^ (where M is 

the mass of a projectile), but is rather given by (3c, where 
/32 = 1 - 1/^2 = 1 _ [Mc^/{Mc^ + r))2, and 7 is the 
Lorentz factor of the projectile. 

Furthermore, Fi and F2 in ([3]) are given by 



Fi{v) = A, 
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The corresponding fitting parameters taken from Ref . |15j , 
Al ... El, A2 ... I?2, ct are listed in tables I and II. The 
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Table 1. Fitting parameters for the three outer shells of water 
molecule with the ionization potentials 7i = 12.61 eV, I2 = 
14.73 eV, h = 18.55 eV. 



^1 


Bi 


Ci 


Di 


El 


A2 


B2 


C2 


D2 


a 


0.97 


82 


0.4 


-0.3 


0.38 


1.04 


17.3 


0.76 


0.04 


0.64 



Table 2. Fitting parameters for two inner shells of water 
molecule with the ionization potentials I4, = 32.2 eV, /s = 
539.7 eV 
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Bi 


Ci 


Dr 


£1 


A2 


B2 


C2 


D2 
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1.25 


0.5 


1.0 


1.0 


3.0 


1.1 


1.3 


1.0 


0.0 


0.66 



cut-ofF energy u)"^^^ is given by 



wf ^'^ = 4z;2 . 






(6) 



where the first term on the right hand side represents the 
free-electron limit, the second term represents a correction 
due to electron binding, and the third term gives the cor- 
rect dependence of the SDCS for w < 1 [15]. For w » 1, 
Rudd's formula would have the correct relativistic asymp- 
totics if Fi, given by Q, is replaced by the following ex- 
pression, 



Fi{v)=Ai 
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whose asymptotic behaviour is clearly similar to the well 
known Bethe-Bloch formula for energy loss: 
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where (/) is the average ionization potential for water 
molecule. This correction reveals itself as an increase of 
the cross section at high energies. In our calculations, we 
omitted the term —0^ in the numerator of Eq. ([7]) after 
checking its negligible contribution in the energy range of 
interest. The SDCS for water as a function of the projec- 
tile's and secondary electron's energies is plotted in Fig.[TJ 
In our previous consideration [13] . the charge of the 
projectile, z, which enters into the expression ([3]) was 
taken to be constant and equal to the charge of a fully 
stripped projectile nucleus. This led to a prediction of un- 
physically large heights of the Bragg peak in disagreement 
with experimental data. This indicates that the charge 
transfer, or electron capture by a charged projectile, takes 
place in the vicinity of the Bragg peak, effectively reduc- 
ing the charge of the projectile [17]. Therefore, z in ([3]) 
should be replaced by an effective charge Zeff which de- 
creases with decreasing energy making the cross section 
effectively smaller. In this work, we use z^ff given by 
Barkas [IH] which represents the average ion charge at 
a given ion velocity j3c: 



z^ff = z(l - exp(-125/?z-2/3)) 



(9) 
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Fig. 2. Average energy of secondary electrons produced as 
the result of impact ionization as a function of kinetic energy 
of ^^0^+ ions. 



The effective charge given by this expression slowly changes 
at high projectile velocity, but rapidly decreases in the 
vicinity of the Bragg peak. As a result, the charge trans- 
fer does not affect the position of a Bragg peak, but it 
may significantly affect its height. 

Integration of the singly-differentiated cross section over 
secondary electron energy W gives the total cross section 
of impact ionization by the ion with the kinetic energy T: 



.m^/ '-^^n. 



(10) 



This quantity is important for all our calculations. Infinity 
as the upper limit of the integration in this and other ex- 
pressions that follow is not physical, but still appropriate 
because of the exponential factor in denominator of ([3]), 
which effectively truncates the integral, [15] . 

The next characteristic that we can obtain from the 
SDCS is the average energy of the secondary electrons, 
which is given by 



W^ave(T) 
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The result of this integration is shown in Fig.dj It is worth 
noting that on the average, the secondary electrons are 
not very energetic, about 50-60 eV, and the average en- 
ergy levels out as the energy of the projectile increases. 
This fact is important for estimating the probability of 
avalanches initiated by the first generation of emitted elec- 
trons. 

The stopping cross section, defined as 



O-st = 



pOQ 



u^^W-^- <-) 



gives the average energy lost by a projectile in a single 
collision, which can be further translated into energy loss 
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Fig. 1. Singly differentiated ionization cross section for ions 12q6+ interacting with water as a function of the projectile's 
kinetic energy T and the energy of ejected electrons W . 



400 


1 1 1 


1 1 1 1 


300 


- ^"\^^ 


_ 


3 
> 

S 200 




\ : 


H 




\. 


100 


- 


\ - 




Ill 


1 , 1 i\ 



50 



100 



150 
X (mm) 



200 



250 



300 



Fig. 3. Energy of ^^C^"'' ions as a function of the penetration 
depth in water. The initial ion energy is 400 MeV per nucleon. 



within an ion's trajectory segment, dx: 

-nast{T) . 



dT_ 
dx 



(13) 



This quantity is known as the linear energy transfer (LET), 
one of the most important quantities in radiation biology. 
The LET found from Eq. dni) is a function of the ki- 
netic energy of the ion rather than the ion's position along 



the path in the medium. The dependence of LET (and 
other quantities) on this position, however, is more suit- 
able for cancer therapy applications. Integrating inverse 
LET, given by ((T3|), we get 



KT) = 



dT' 



\dT'/dx\ 



(14) 



where Tq is the initial energy of the projectile. We obtain 
the correspondence between the position of the ion along 
the path and its energy. The dependence x{T) for carbon 
ions is shown in Fig. [31 

This figure helps to remap all quantities of interest so 
that they depend on x rather than on T . The depth de- 
pendence of the average LET as a function of x is shown 
in Fig. m In this figure, our calculated LET is compared 
with the predictions of the MCHIT model [8|9|10j based 
on the GEANT4 toolkit [11I12J . The Monte Carlo cal- 
culations were performed both with and without taking 
into account nuclear fragmentation reactions. Experimen- 
tal data by Schardt et al. il6j are also shown for compar- 
ison. Our calculated position of the Bragg peak is about 
5 mm (2%) deeper than predicted by simulations and ob- 
served by the experiment. This discrepancy arises from 
neglecting the projectile energy loss due to non-ionization 
processes when the energy is spent for the excitation of 
neutral atoms without ejection of electrons. Indeed, it is 
estimated [TJ as 5-15% of the total projectile energy loss. 
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Fig. 4. Linear energy transfer of 400 MeV/u ^^C^''' ions as a function of penetration depth in water. The MCHIT resuhs with 
and without nuclear fragmentation are shown for comparison. Points - experimental data [16J 



So far our consideration neglects the fact that the pro- 
jectile energy loss is a stochastic process, and \dT/dx\ was 
calculated as the average LET. Therefore, we introduced 
the energy-loss straggling which essentially improves the 
agreement with the MCHIT results obtained taking into 
account the stochastics of energy-loss process. In calculat- 
ing the energy-loss straggling, we followed the prescription 
of a semi-analytical model described in Ref. [21] , 
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where xq is a maximum penetration depth of the pro- 
jectile, and astr — 0.8 mm is the longitudinal-straggling 
standard deviation computed by HoUmark et al. |22j . for 
a carbon ion of that range of energy. The results are pre- 
sented in Fig. [5l In order to facilitate the comparison of 
our analytical model with MCHIT, in Fig. [H MCHIT re- 
sults are presented without taking into account nuclear 
fragmentation reactions, and our data are plotted with a 
proper shift in x, accounting for the difference between 
the peak positions. This allows us to consider separately 
the effects of electromagnetic interactions of beam nuclei. 



As confirmed by MCHIT simulations [12] , nuclear frag- 
mentation reactions become important for heavy-nuclei 
beams and deeply-located tumors. For example, both ex- 
perimental data 16] and MCHIT calculations |20] indicate 
that more than 40% of primary 200 MeV/u 12^6+ nuclei 
undergo fragmentation before they reach the Bragg peak 
position, and this fraction exceeds 70% for 400 MeV/u 
12C6+ beam. 

What are the effects of nuclear reactions? First, the 
beam is attenuated since ions are leaving the beam. Sec- 
ond, new projectiles are formed and they may have dif- 
ferent penetration depths than the original carbon ions. 
Most important, protons and a-particles, which dominate 
the secondary products of such reactions [IS] have much 
longer penetration depths at the same energy (per nu- 
cleon). This results in a substantial tail after the Bragg 
peak shown in Fig. |4l 

The next quantity of interest is the number of sec- 
ondary electrons with kinetic energy VF, produced by a 
single ion on the interval Ax at the depth x correspond- 
ing to a certain kinetic energy T of the ion. This quantity 
is a product of the remapping of Eq. ([T]) integrated over 
solid angle; it can be written as follows: 
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Fig. 5. Linear energy transfer of 400 MeV/u C ^ ions as a function of penetration deptli in water with and without straggling 
compared with the MCHIT results without nuclear fragmentation. (*)Our results are shifted by 4.6 mm for convenience in 
comparing the peak shapes. 
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Fig. 6. Number of secondary electrons with kinetic energy W , 
produced by a single ion in a 1-nm interval vs. kinetic energy of 
the secondary electrons, for different ion's energies T(MeV/u) 



Four curves shown in Fig. [6] show this dependence at dif- 
ferent energies or depths of the ion. This number, however, 
represents only the secondary electrons produced by the 
projectile. In order to determine the total number of elec- 



trons in the avalanche, one has to analyze the ionizing 
capabilities of these secondary electrons. 

To address this issue, we calculate the fraction of sec- 
ondary electrons with energies higher than the threshold 
of ionization produced in a single collision. 



1 



a{T) 



dW' 
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The dependencies for different electronic shells are shown 
in Fig. [T] These results are somewhat surprising since they 
indicate that the fraction of secondary electrons capable of 
further ionization decreases with increasing energy of the 
ion. This is probably a result of averaging; even though 
higher-energetic secondary electrons are produced, there 
are many more lower-energy ones (see Fig.[T|). The fraction 
of secondary electrons capable of further ionization is sub- 
stantial for secondary electrons produced by the projec- 
tile, even at relatively low projectile energies; however, sec- 
ondary electrons produced by secondary electrons are even 
less energetic, so that further ionizations are much less 
probable. This is supported by calculations using Rudd's 
model for incident electrons. However, it can also be ex- 
plained by a simple estimate: Fig. [2] indicates that the av- 
erage energy of secondary electrons in the vicinity of the 
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Fig. 7. Fraction of secondary electrons produced by ^^C^+ 
capable of further ionization. Different lines indicate different 
electronic shells of the water molecule. 



Bragg peak is about 45 eV. The maximum average energy 
that can be transferred to the next generation secondary 
electron is just (45 — Ii)/2, which is just over 12 eV for 
the outermost electron, an energy barely enough to cause 
further ionization. Therefore, we conclude that there is 
no avalanche of secondary electrons and that Eq. (|16p 
gives a reasonable estimate of the number of secondary 
electrons, within a factor of 2. Thus, we have analyzed 
the track of a single ion: we obtained the LET (depicting 
the Bragg peak, its shape and position), we obtained the 
energy spectrum of the secondary electrons, and we esti- 
mated the avalanche effects related with ionization caused 
by secondary electrons. 

3 Direct consequences of ionization and DNA 
damage 

In this section we estimate three components that con- 
tribute to the DNA damage: an electron plasma, free rad- 
icals produced as a result of ionization of water molecules, 
and heating of the medium caused by ion passage. 



3.1 Electron plasma and free radicals 

Let us estimate the number density of the electron plasma 
produced by an ion. In Section [21 we found the secondary 
electron energy spectrum and estimated the effect of fur- 
ther ionization. In order to deduce the density of the elec- 
tron plasma, we have to estimate the volume occupied by 
these electrons. The crucial quantity for this estimate is 
the penetration length of the secondary electrons caused 
by the projectile. In principle, this estimate can be done in 
the same fashion as it was done for the ions in the previous 
section; however, the energies of the secondary electrons 



are such that Rudd's analysis may not be reliable |15j . 
Therefore, we will try to rely on other available sources. 

Direct experimental observations of penetration depths 
of electrons in liquid water are scarce. There are values 
recommended by the International Commission on Radia- 
tion Units and Measurements for 10 - 1000 keV electrons 
23 24J, and there are data for thermalization lengths of 
"subexcitation electrons" for lower energies, 0-4 eV |25j . 
However, there is a gap in experimental data for electron 
energies, which are the most interesting from the view- 
point of cancer therapy applications. 

The average energy of secondary electrons in the Bragg 
peak area (T ~ 0.3 MeV per nucleon) is about 45 eV 
(see Fig. [2). Even though there are no direct experimen- 
tal values for the penetration lengths of electrons in this 
energy range, interpolation of existing experimental data 
points [55] and theoretical predictions based on Monte- 
Carlo simulations [27] suggest that the penetration lengths 
of such electrons are of the order of 10 nm or less. Com- 
parison of these lengths and typical distances between the 
ions in the beam leads us to conclude that under realistic 
conditions ion tracks do not interact. That is, every ion 
acts independently and the total biological effect equals 
the sum of the damages from individual ions, and therefore 
the collective action of ions in the beam can be neglected. 

In order to get realistic order of magnitude estimates, 
we assume that the energy deposited by an ion is thermal- 
ized within a tube of a radius r ^ 10 nm. 

The total ionization cross section for carbon ions i2q6+ 
of kinetic energy T ^ 0.3 MeV per nucleon is about a « 
40 • 10~^^ cm^ = 0.4 nm^. The majority of secondary elec- 
trons are produced in the following process: 



,(.b 



H2O 



C 



6+ 



H20^ 



(18) 



The number of electrons produced per one nm of the ion's 
trajectory is then na = 33 nm~'^ x 0.4 nm^ w 13 nm^-'^. 
This quantity is easily translated into the number density 
of the electron plasma. 



13 nm 
ttIOO nm^ 



0.045 nm~3 = 0.45 • 10^" cm"^ 

(19) 



Both secondary electrons and ionized water molecules 
react to form free radicals 



H20^ 



H^ 



H2O 



H 



OH' 



OH 



(20) 
(21) 



This means that the induced concentrations of free radi- 
cals can be of the same order as the concentration of the 
electron plasma. 

OH radicals play an important role in damaging DNA. 
They interact with nucleobases damaging them; moreover, 
"secondary" free radicals may result from these interac- 
tions ^. Further analysis is needed in order to obtain 
quantitative estimates of the yields of DSB's and other 
types of irreparable DNA damage due to these effects. 
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3.2 Dissociative recombination by low energy electrons 

Traditionally, it had been believed that non-thermalized 
secondary electrons can not cause any significant geno- 
toxic damage, unless they are already solvated in water 
and can participate in chemical reactions of the type (j2ip 
or unless they have enough kinetic energy. The singly dif- 
ferentiated cross section, however, falls rapidly with en- 
ergy of the ejected electron and therefore energetic sec- 
ondary electrons are not abundant, cf. Figs. [T] and [2l 

In 2000, Sanche et al. [28] argued that 1- to 20-eV 
secondary electrons "induce substantial yields of single- 
and double-strand breaks in DNA, which are caused by 
rapid decays of transient molecular resonances localized on 
the DNA's basic components" [28] . In their experiments, a 
monoenergetic electron beam at various incident energies 
was used to irradiate solid DNA samples. They observed 
that one electron out of 1000 - 2000 (depending on the 
electron energy) induced a single strand break and one 
electron out of 5000 - 10000 induced a DSB of a DNA 
molecule [29j . 

This type of experiment, however, can not be directly 
related to the problem of DNA damage induced by ions. 
The difference is in the density of free electrons. In Sanche's 
experiments the electron density was of the order of 10** cni^ 
16 orders of magnitude lower than is created by an ion in 
the Bragg's peak area, as it is estimated in the previous 
section. 

One can nevertheless estimate the average number of 
electrons created in a volume occupied by one convolu- 
tion of the DNA molecule. The DNA radius is about 1 
nm, while the convolution (linear) length is about 4 nm, 
so that the volume of interest is about Vdna = lO^^*' cm"^. 
Assuming that at least two electrons are needed in the vol- 
ume of one DNA convolution to produce a double strand 
break, one can estimate that the required electron density 



in the Bragg peak region), the tube's radius r and the 
density of water p = 1 g/cmT^: 



^DNA 



2 • lO^'' cm" 



(22) 



It is remarkable that such a density is of the same order 
of magnitude of the one that has been predicted by our 
simple estimates in the previous section. 



3.3 Local heating 



Q = pur AxcAT, 



AT ■ 



Q/Ax 
pnr'^c 



(24) 
(25) 



Substituting the numerical values, one obtains AT ^ lOK 
for r = 10 nm and for the lower estimate of the tube's 
radius, r = 3 nm, we find AT '^ lOOisT. 

This estimate shows that the local heating within the 
electron thermalization radius along the ion track can 
be quite substantial. The melting temperature of a DNA 
molecule is about 85° C. Therefore, in our opinion, there 
is a need for a more thorough study of the local heat- 
ing mechanism which would include investigation of such 
problems as the detailed description of the thermalization 
of secondary electrons in water, heat transfer in water, 
and modeling of DNA dynamics under local heating. To 
the best of our knowledge, this local heating mechanism 
has not yet received any attention. 



4 Conclusions 

We have presented an approach to modeling ion-beam 
therapy having considered the effects initiated by an ener- 
getic projectile such as a carbon ion incident on biological 
tissue (represented by liquid water) . We analyzed passage 
of the ion through this medium taking into account the 
main processes that cause energy loss by the projectile. 
We succeeded in making quantitative predictions of effects 
such as the energy spectrum and abundance of secondary 
electrons and local heating caused by the projectile. Then 
we built our estimates of DNA damage on these predic- 
tions. In principle, the final estimates may be related to 
the energy deposition by the projectiles, but the micro- 
scopic analysis of the whole scenario is vitally important 
in order to provide more accurate predictions that may 
eventually contribute to the protocol of cancer therapy. 
Then they can be formulated into the language of dosages, 
energies, radiation rates, etc. Before that, more research 
is needed in understanding the mechanisms of DNA dam- 
age on a microscopic level and relating these mechanisms 
to the physical parameters of secondary electrons and free 
radicals produced by the ion beam. 



Finally, let us estimate the local heating produced in a 
tube of 10 nm radius by a carbon ion of 0.3 MeV per 
nucleon. The simplest way to do that is to use the ther- 
modynamic equation 



Q = pcAT 



(23) 



relating the heat transferred, Q, to a system with the sys- 
tem's mass p, specific heat capacity c (4.2 J/g K for wa- 
ter) and the increase of system's temperature AT. AT can 
then be found from the linear energy transfer Q/Ax (typ- 
ical value for ^^C ions is 100 MeV/mm = 1.6 lO'^'^J/cm 
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